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Abstract 

We present a comprehensive study of the phase transitions in the single-field 
reaction-diffusion stochastic systems with field-dependent mobility of a power-low 
form and the internal fluctuations. Using variational principles and mean- field the- 
ory it was shown that the noise can sustain spatial patterns and leads to disordering 
phase transitions. We have shown that the phase transitions can be of critical or 
non-critical character. 

Key words: Nonlinear diffusion; Internal noise; Phase transition 
PACS: 47.54.-r, 05.40.-a, 05.70.-a 



1 Introduction 



The ability of noise to induce spatial ordering has received a special attention 
in last two decades. It is well known that in many situations noise can actually 
play a constructive role, to name just a few one can consider: noise- induced 
transitions in zero-dimensional systems [1] , noise induced phase transitions in 
extended systems [2], patterns formation [3], coupled Brownian motors [4,5], 
etc. Most of the noise-induced effects are caused by an external fluctuations 
and attributed to a short time instability. In particular, noise-induced patterns 
and the phase transitions usually have a dynamical origin. It is principally im- 
portant that in the case of the internal noise, obeying fluctuation-dissipation 
relation, a spatio-temporal coherence in nonlinear systems can not be observed 
in a short-time scale. Such ordering processes follow an entropy-driven mech- 
anism when a kinetic coefficient /mobility is a function of a stochastic field. It 
leads to the fact that a stationary distribution is described by a free energy 
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functional reduced to a Lyapunov functional for deterministic dynamics and 
by an entropy contribution related to the field-dependent mobility [6,7,8,9]. It 
was shown by the most of works concerning such phenomena that quantitative 
change of the system behaviour is caused by the entropy variability. 

It was shown that internal fluctuations can sustain the spatial coherence (pat- 
terns) in nonlinear react ion- diffusion stochastic systems [8,10]. In the simplest 
case such models have two essential features: the local dynamics and a trans- 
port phenomena. First of them is determined by the chemical reactions in the 
system and the second one relates to a diffusion of chemical species. Introduc- 
ing a fluctuating source into the model, we arrive at the stochastic descrip- 
tion of non-equilibrium system behaviour with possible spatial order. In this 
work we consider a prototype model of the reaction-diffusion systems where 
an internal fluctuations are of a multiplicative character (the noise intensity 
depends on a concentration field) (see Ref.[8] and citations therein). In the 
most of works devoting to study the system behaviour with a field-dependent 
mobility, the corresponding functional dependence is as follows: in the mixed 
state the kinetic coefficient is large, whereas in dense ones it has small val- 
ues. A deterministic analysis of the systems with a generalized field-dependent 
mobility was performed in Ref. [11], where the form of the kinetic coefficient 
was proposed by authors solely. In our consideration we use a generalized 
approach basing on a nonlinear master equation proposed by G. Kaniadakis 
[12] and derive the generalized form for the kinetic coefficient. Introducing the 
corresponding Fickian diffusion term into the model and a term to describe 
a local dynamics of the field variables, we arrive at a deterministic model for 
the reaction-diffusion systems. Considering a system in real conditions we take 
into account a fluctuation source related to an obtained mobility. 

According to the well known approaches of entropy-driven phase transitions 
[6,7] one can investigate a possibility of patterns formation and corresponding 
phase transitions. Our goal in this article is to perform a detailed study of 
spatial system's coherent behaviour with a generalized kinetic coefficient in a 
bistable stochastic system. With help of variational principles we investigate a 
possibility of noise to sustain stationary spatial structures. A picture of noise- 
induced phase transitions will be studied with the mean-field approach and 
computer simulations. 

Our paper is organized in the following manner. In the next section we present 
the general formalism to investigate the noise-induced phase transitions and 
patterns formation. In Section 3 we derive an expressions for the field-dependent 
diffusion coefficient using the formalism of nonlinear master equation and de- 
formed Boltzmann-Gibbs statistics. According to the presented formalism the 
noise-sustained structures are considered in Section 4. The noise-induced phase 
transitions with the entropy-driven mechanism are studied in Section 5. Main 
results and prospects for the future are presented in the Conclusions. 
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2 General theory 



Our starting point is the evolution equation for a particle density p = p(r, t) 
written in a most general form of the continuity equation [12] 



dp 
di 



Dl (p) d -^Wp 



(1) 



Here we take into account diffusion flux only with concentration-dependent 
effective diffusion coefficient assumed in the form 

DM^D^p)^^-, D = const. (2) 

From a formal viewpoint the nonlinear diffusion equation (1) can be de- 
rived from a master equation assuming that probabilities of the microscopic 
transitions explicitly depend on the concentrations of the initial and arrival 
states, described by functions fi(p) and z^(j»), respectively. It was shown (see 
Refs.[13,14]) that functions 7(jo) and n(p) in Eq.(l) are defined as follows: 
lip) = f i ip) u ip)y K iv) — MpVKp)- A 11 explicit form for pairs (/i, v) or (7, ft) 
can be set due to the physics of the problem under consideration. 

Introducing chemical reactions described by a rate f(p), a generalized conti- 
nuity equation for the concentration of particles reads: 

^ = /(p) + V [£> e /(p)Vp] . (3) 

Formally, Eq.(3) can be written with help of the Lyapunov functional for the 
deterministic dynamics. Indeed, introducing notation for such a functional in 
the form 

F = J dr {±[D ef (p)Vp] 2 + <p(p)} (4) 

with 

v 

<P=-J fip')D ef (p')dp', (5) 


the deterministic evolution equation (3) takes a variational form [8] 

dp = l_5_F\pl 

dt D ef (p) Sp ' 1 ] 



Considering the system in realistic conditions one needs to introduce a fluctu- 
ation source £(p;r, t). In our stochastic analysis we assume that such fluctu- 
ating source obeys the fluctuation-dissipation relation [15] and has following 
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properties: 



2(T 2 

<£(p; r, t)) = 0, (£(p; r, f)£(p; r', f)> = 7m ^ " " ^ ^ 

where a 2 is an intensity of the corresponding internal multiplicative noise. For- 
mally, one can introduce an external noise related to fluctuations of a control 
parameter addressed to a local dynamics. In this paper we study influence of 
an internal fluctuations source only on the phase transitions picture with the 
concentration-dependent diffusion coefficient and the chemical reactions. In 
the following analysis we use the Stratonovich interpretation of the Langevin 
equation 

ap _ 1 ,F H + 1 



dt D ef (p) 5 P yjD ef {p) 

where £(p;r,t) = [D ef {p)]- l / 2 l(v,t). 

Considering stationary properties of the system we exploit a stationary prob- 
ability density functional, obtained as a solution of the corresponding Fokker- 
Planck equation [16]. In the framework of standard technique such stationary 
functional takes the form [2,6,17] 

P st oc exp 
the effective energy functional 

U ef \p]=F\p]-£fdx\nD ef (p) (10) 

is defined through the free energy functional F\p\ and the effective diffusion 
coefficient D e f(p). Due to the stationary probability density functional (9) has 
an exact form, the effective potential (10) can be used to study possibility of 
the structure formation under the multiplicative noise influence. Moreover, 
we can apply the mean-field theory formalism to investigate the noise-induced 
phase transitions in systems of such a kind. 



3 Model 



The first problem we deal with is to set a generalized form for the effec- 
tive diffusion coefficient D e f(p), defined according to Eq.(2). The main crite- 
rion for the function D e f(p) is a bell-shaped form in the interval p G [0, 1]. 
We assume that following properties are satisfied: _D e /(0) = D e f(l) = 0; 
D e f (1/2) = D^ ax \ where £)(™ ax *> [ s a maximal value. It means that fluc- 
tuations are large in a mixed state characterized by the value p = 1/2; in 
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dense states (p — 1, p — 0) no fluctuations are realized. Considering a gen- 
eral problem, related to a complex systems investigation, let us assume that 
the function D e f(p) has a power-law form. Indeed, as it was shown before 
the complex systems can be described by the nonlinear continuity equation 
and are characterized by a power-law form for main statistical quantities (see 
Refs. [18, 19,20,21]). Most of complex statistical systems are described by de- 
formed Boltzmann-Gibbs statistics or q-statistics exploiting deformed loga- 
rithm's and exponential functions [21]. By now there are two well known kinds 
of deformations, proposing by G.Kaniadakis and C. Tsallis (see Refs. [13,21], 
respectively). In this paper we operate with a mathematical construction of 
the q-deformed logarithm lnx — > \ia q x = {x 1 ~ q — 1)/(1 — q) [21], the expo- 
nent q is a non-additivity parameter; in the limit case q = 1 one arrives at 
usual logarithm and the standard Boltzmann-Gibbs statistics. Thus, inserting 
q-deformed logarithm into Eq.(2), one has 

D ef (p) = £^(p)K(p)-«^. (11) 

In the most of works related to the nonlinear continuity equation (1) the 
functions 'yip) and nip) are assumed to be linear versus its argument. In the 
simplest case one has: k(p) = p, j(p) = p. It results to the construction of the 
form 

D ef (p) = Dp 1 '*. (12) 

It is seen that the linear approximation of j(p) and n(p) does not give the 
bell-shaped form for the effective diffusion coefficient. Evolution of the system 
with the effective diffusion coefficient (12) was described in Ref. [22]. It was 
shown that such complex systems manifest an anomalous diffusion. 

To satisfy the main condition for the function D e f(p) let us assume nonlinear 
constructions for 'yip) and n(p). According to the fact that complex systems 
are self-similar, let us write 'yip) in the form 

jip)=p a , (13) 

where the exponent a G (0, 1). To define the function we suppose that the 
derivative d«(p)/dp in Eq.(ll) gives a power-low dependence, i.e. d«(p)/dp = 
k@, j3 > 0. Such an assumption determines the power- low form for the function 
k(p). Indeed, after some algebra one can obtain k{p) — [C(/3 — 1) + (1 — j3)p}~ , 
where C > is a constant. Next, introducing a positive constant C\ = C(/3— 1) 
with f3 > 1 one can represent the function n(p) in the form of the Tsallis ex- 
ponent [23]: 

«(p) = [d + (1 - (3)p\^ = ex P/3 (14) 

The exponent j3 plays a role of the non-additivity parameter. The Tsallis 
exponent in Eq.(14) becomes the usual one at (5 — 1. Substituting expressions 
for the functions 7(p) and from Eq.(13) and Eq.(14), respectively, into the 
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equation (11) and performing trivial calculations one arrives at the power-low 
form for the effective diffusion coefficient: 



D ef (p)=Dp a [C 1 -((3-l)p} 5 , 5 



q-P 
0-1' 



(15) 



The function D e f(p) has a bell-shaped form only if 5 > or if q > (3 with f3 > 1. 
Formally, the effective diffusion coefficient D e f(p) defined by the exponent 
a 7^ 5 is a non-symmetrical function with respect to the point p = 1/2. In 
further investigation we are interested in studying the system properties when 
D e f(p) is the symmetrical function. To this end we assume 5 — a. It allows 
to obtain the relation between the both non-additivity parameters q and (5 in 
the form 

(5<q<2(3-l. (16) 

Therefore, we arrive at the symmetrical form for the effective diffusion coeffi- 
cient 

p 



Def( P ) 



A) a 



1 - 



Ps 



(17) 



where p s = C\/{f3 — 1) is a saturation concentration, (3 = 2DC" is a constant. 
The obtained formula (17) can be derived directly if probabilities fx(p) and v(jp) 
are known initially. Indeed, using the approach developed in Ref. [14] one can 
express functions 'yip) and nip) through fi(p) and v(p), and after find the 
alternative construction for the effective diffusion coefficient: 



Defip) = D 



Ap). 



1-9 r 



d d 
v—H-H—v 
dp dp 



As it follows from our consideration densities /i(p) and u(p) can be obtained 
using relations between j(p) and n{p). After a trivial algebra one find: 



Mp) = b(p) K (p)] 1/2 , v{p) = [i(p)/k(p)} 1/2 . 



(19) 



Therefore, the form for the effective diffusion coefficient is well defined. 



Next, to derive a model for the function f(p) that describes possible chemical 
reactions in the system. Let us assume that in the deterministic regime there 
are three stationary concentration values p^\ p^\ p^ with p^ ^ p$ \ i ^ j. 
Thus, using the theory of dynamical systems one can suppose the deterministic 
force in the form f(p) = — Oi(p — Pa?)- I 11 the simplest case the following 
construction can be used: 



f(p)=e(p-\)-(p-\) 3 , (20) 

where e is a control parameter, A determines the equilibrium concentration 
magnitude. In such a case stationary states are determined by the values: 
Pq 1 -* = A, Po 2,3 -* = A ± yfi. In the following analysis we choose A = 1/2. 
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Next, let us show that the corresponding internal multiplicative noise leads 
to a short-time instability of the mixed/disordered state p — 1/2. The linear 
stability analysis can be performed for an auxiliary field y(r,t) = p(r,t) — A. 
The corresponding linearized Langevin equation takes the form 

= £ y + ^y + ?y + ^(r, t) , (21) 

where x 2 = l/2p \ 2a , a 2 = (2aK7 2 )/(A)A 2(a+1) ). In the Fourier space Eq.(21) 



is 



dV{± ^ t] = ^(± k > *) " * 2 k 2 y(±K t) + a 2 y(±k, t) + x-^(±k, t). (22) 
Writing the dynamical equation for the two point correlation function 

^<y(M)y(-M)> = ^(k,i)J^(-k,^ + (y(-M)^y(M)^ (23) 

and calculating correlators (£-kyk), (6<y-k) with the help of Novikov's theorem 
[24], one arrives at a dynamical equation for the structure function S(k,t) = 
(y k (0 y-k(O) in the form 

o 2 

-S(k, t) = 2(e - x 2 k 2 + a 2 )S(k, t) + 2—. (24) 
at x 

It is seen that the internal multiplicative noise leads to a short-time instability 
of the mixed state p — 1/2. The stationary value of the structure function is 



X' 



Hence, the homogeneous state p(r,t) = 1/2 is stable only if e < >c 2 k 2 — a 2 



4 Noise sustained structures 



As was shown in Ref. [8] the extrema of the effective potential U e f\p] corre- 
spond to the stationary noise-sustained structures p st . Such structures can 
be computed as an solution of equation 5U e f\p\/5p = 0, where the effective 
functional (10) has the form 



Ueflp] = J 



dr 



<p(p) + \ {D ef {p)Vpf 



-^|drln/J e/ (p). (26) 
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Making the first variation of U e f[p] with respect to p equal to zero, we arrive 
at the equation for the stationary structures 



Ap = - 



a 2 dD ef (p) t f(p) 



+ 



+ 



2D 3 ef (p) dp D ef (p) D ef (p) dp 



9Def(p) (Vp)2 



(27) 



P=Pst 



Considering homogeneous states, one can put Vp = Ap = in Eq.(27). The 
reduced equation gives the most probable stationary states p e , defining ex- 
trema positions of the function U e f(p). The corresponding solutions of such 
equation are shown in Fig. la. It is seen that with an increase in the noise 
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Fig. 1. Most probable values of the quantity p versus the noise intensity a 2 at e = 1.0 
(a) and phase diagram of a system at (3$ = 2.0 (b) 

intensity two different values of the most probable concentration degenerate 
at a bifurcation point a 2 = of. Despite the fixed point p e = 1/2 exists al- 
ways, the two additional fixed points (upper and lower curves in Fig. la) are 
observed at a 2 < a 2 . At a 2 > a\ the stationary probability function has one 
extremum only located in p e = 1/2. The phase diagram a (of) illustrating the 
dependence of the bifurcation point position at various values of e is shown in 
Fig. lb. One can see that with an increase in the exponent a the bifurcation 
value <7 2 decreases. Therefore, if D e f decreases sharply at p ~ and p ~ 1, 
then the stationary distribution becomes unimodal at large values of the noise 
intensity a 2 . 

Next let us investigate stationary structures p st . Firstly, we consider the vicin- 
ity of the homogeneous solution p e — 1/2. In a linear approximation one can 
put (Vp) 2 = in Eq.(27) and expand the left part of Eq.(27). Then, we obtain 
the second order differential equation in the form 



|§ ~ Akf (p - p st ) , 



A = +1 if a 2 > a 2 , 
A = -1 if a 2 < al 



(28) 



where h = ^{e-a 2 )/x 2 . So, one can conclude, that at a 2 < a 2 the station- 
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ary state p e = 1/2 is locally stable and represents a center in the phase space 
(p, Vp); in the opposite case (a 2 > of) it changes the stability and becomes 
a saddle. It is seen that the dependence ki{a 2 ) has monotonically decreas- 
ing character. As it follows from our consideration the stationary periodic 
structures are formed in the vicinity of the point related to maximum of the 
function U e f(p). 

In the numerical studying of the noise-induced spatial patterns we integrate 
Eq.(8) on a rf-dimensional lattice of the mesh size /. In the discrete space 
Eq.(8) can be written as follows: 

jr'grm+(j) " l j£nn(i) 



v 2 dDe f (Pi) 1 + 1 c 
2 d Pl D 2 ef ( Pi ) fjjfcj ' 



(29) 



where % — 1, . . . , N d enumerates the element of the square lattice of N d cells; 
periodic boundary conditions are used. The second term in Eq.(29) represents 
approximation of the gradient |Vp| 2 , where nn + (i) indicates nearest neighbors 
in the positive direction of each axis, whereas the third term is related to a 
discrete Laplacian (Vp — > J2 ^ijPjj nn{i) denotes nearest neighbors of 

j£nn(i) 

the site i). To describe spatial patterns we use a spherical averaging of the 
structure function S{k,t) = J r S(k,t)dT, where T k is a spherical shell of a 
radius k. A convenient formula is 

S(k,t) = ± £ S(k,t). (30) 

JVfc k<k<k+Ak 

All calculations were performed in a two-dimension square lattice of 120 x 120 
cells with lattice scaling / = 1, and integration time step r = 5 x 10 -3 . The 
stationary spherically averaged structure function is shown in Fig. 2 at different 
values of the noise intensity. 

It is principle, that an increase in the noise intensity a 2 results to a shift of 
the peak of S(k) to the small values of the wave vector k. This is different 
from what happens at the early stages of evolution. This fact indicates that 
when the system starts to evolve a multiplicative noise leads to instability 
of the state p — 1/2; the corresponding values of the wave vector k increase 
with the noise intensity growth. For large times (in the stationary case) the 
situation is quite different. Here with an increase in the noise intensity the 
peak of S(k) moves toward smaller values of k. This fact is related to an 
entropy mechanism of the phase transitions. Indeed, the effective energy func- 
tional (10) is defined through the initially known free energy functional F[p] 
and a contribution — a 2 /2 J dr In D e f(p). The last one gives an effective en- 
tropy S e f = —1/2 / dr\n D e f(p) multiplied by the noise intensity a 2 . In such a 
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Fig. 2. Dependence of structure function S vs. wave vector k at e = 1.0 a = 0.5, 
Po = 2.0 and various values of noise intensity 

case we arrive at the standard thermodynamic definition of an internal energy 
U e f = F + & 2 S e f. It is well known that in the theory of the entropy-driven 
phase transitions self-organization processes are not related to the short-time 
instability of the mixed/disordered state [6], its caused by entropy variations 
following from concentration-dependent mobility. In our case we have a quite 
similar situation. In particular, at early stages of the system evolution a noise 
destabilizes the disordered state p — 1/2, whereas at t — > oo the entropy con- 
tribution plays a crucial role: it leads to patterning with small values of k. 
Obtained numerical results are in good corresponding with analytical predic- 
tions: in the stationary case the peak of S(k) moves toward small values of 
the wave vector when the noise intensity increases. 



5 Mean Field Theory 

In this section we will consider the noise-induced phase transitions with the 
above entropy mechanism. To this end we construct the Weiss' mean-field 
(MF) approximation, based on the stationary distribution function. Let us 
approximate the gradient term in the functional (10) by the sum over nearest 
neighbors on the lattice as Eq.(29) shows. In the MF approximation we replace 
the exact value of the neighbors by a mean field rj = (p). It leads to the relation 
(Vp) 2 — > (p — (p)) 2 = (p — i]) 2 . In such a case the quantity rj can be used as 
an order parameter. In this procedure we neglect fluctuations in neighboring 
sites. The value of the order parameter can be computed self-consistently as 
follows: 



where P s t(p',v) i s the stationary distribution in the MF approximation de- 
pending on the mean field 77. As usual, the parameter rj allows to identify a 
phase transition from a disordered state with rj = i] c to an ordered one char- 
acterized by the value rj 7^ i] c ; the symmetry that leads to r] = r\ c is embedded 




(31) 
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in the system. According to the discussion above, the effective potential has 
the form 



1 2 

U ef (p; 77) = -D 2 ef (p)(p -Vf-Y ^D ef (p) " / f{p')D ef (p'W (32) 



and depends on the mean field r\. The stationary probability density function 
is 

P st (p; 77) = Z{rj) exp (-U ef (p; r,)/a 2 ) , (33) 

where Z{j]) takes care of the normalization condition J P st (p;r))dp = 1. Ac- 
cording to formalism proposed in Ref . [7] and quantitative analysis of the func- 
tion $(77) one can expect the second-order phase transitions in a system under 
consideration. Using Newton-Raphson condition, transitions between disor- 
dered and ordered phases occur if 



d$(77) 



drj 



1, 77c = 1/2. 



(34) 



V=Vc 



Introducing the susceptibility in the form 

d$(7/) 
di] 



X 



(35) 



we relate the critical points to the values of a 2 when x — 1- The corresponding 
phase diagram is shown in Fig. 3. 




0.2 0.4 



0.6 0.8 1.0 



Fig. 3. Phase diagram at s = 0.1 (shaded domain corresponds to the non-critical 
phase transitions) 



Here the surface plotting with the help of the solid lines corresponds to the 
critical values of the system parameters, when the second-order phase transi- 
tion are occur. Shaded domain corresponds to the phase transitions which are 
of non-critical character (such situation is observed in magnetic systems when 
a weak external field is introduced into the system). Here the order parameter 
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decreases continuously and has not singularity in the vicinity of the transition 
point related to the shaded domain. The related susceptibility has the broad 
peak in the point corresponding to the system parameters values of the shaded 
domain shown in Fig.3. To understand a principle change in the system be- 
haviour one can estimate contributions of all nonlinearities appeared in the 
effective potential (32). Here one has a competition between the noise char- 
acterized by D e f(p) and the driving force f(p). Indeed, at a ~ 1 as the term 
related to the force f(p) as the term related to ha.D e f(p) are essential and lead 
to phase transition of a critical character. At a <C 1 the main nonlinearity is 
caused by the force f{p) only. Here the system can be approximately described 
by D e f(p) const. At the a- values related to the shaded domain there are 
two above nonlinear terms having weak nonlinearities than the term related to 
Dlf(p — r)) 2 gives. Therefore, at such intermediate values of a there are no non- 
linearities suppressing the order in the system. A size of the domain for such 
transitions Aa = a — a m depends on control parameters fl and e. Let us in- 
troduce an effective order parameter fj — \l/2 — r]\: in the "disordered" phase 
with p — 1/2 one has fj = 0; in the "ordered" phase fj ^ 0. The dependencies 
of the effective order parameter and the susceptibility versus the exponent a of 
the effective diffusion coefficient D e f and noise intensity a 2 are shown in Fig.4. 
From the Fig. 4a one can see that at a > a and a < a m the effective order 




Fig. 4. The character of changing of the order parameter and the susceptibility at 
£ = 0.1, O = 10.0 



parameter fj decreases critically to zero with the noise intensity growth. In the 
domain a G (a m , aio) the dependence fj((T 2 ) has a smooth falling-down char- 
acter. Therefore, such transitions are non-critical. Figure 4b illustrates peak 
of the susceptibility x that, corresponds to the value x(a, a 2 ) = 1 and deter- 
mines the second-order phase transition point. In the interval of the exponent 
a G (a m , «o) the susceptibility x has the broad peak which is characterizes by 
value x — 1 an d relates to non-critical phase transitions. Thus, as it follows 
from our calculations to obtain the critical (second-order) noise-induced phase 
transitions the exponent a should lie in the interval (0, a m ) LK^o? !)■ 
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To verify our MF results we integrate numerically Langevin equation (8) on the 
two-dimensional lattice according to algorithm presented above (see Eq.(29)). 
For the order parameter we use the formula fj — |l/2 — the generalized 
susceptibility measures fluctuations of the field p in the vicinity of the critical 
points, x — (Ap) 2 = (p 2 ) — (p) 2 , averaging over noise realizations, time, and 
the ensemble is taking into account. Solving the discrete Langevin equation 
Eq.(29) numerically on a two-dimension square lattice of 120 x 120 cells with 
mesh size / = 1, and integration time step r = 5x 10~ 3 , we obtain dependencies 
f)(a 2 ) and x((J 2 ) shown in Fig. 5. It is seen, that the order parameter decreases 




0.02 0.03 0.04 0.05 0.01 0.02 0.03 0.04 0.05 0.06 



a) ° 2 b) ° 2 

Fig. 5. The dependence of the order parameter fj and the susceptibility \ vs - noise 
intensity a 2 at e = 0.1, j3o = 3.0 

critically with the noise intensity growth. The generalized susceptibility has a 
well pronounced peak located at the critical transition point. With an increase 
in the exponent a the peak location of the generalized susceptibility x is shifted 
toward small values of the noise intensity a 2 . Our numerical results are in a 
good correspondence with the MF theory. 

Figure 6 illustrates dependence of the order parameter and the generalized 
susceptibility at the system parameters values related to the shaded domain 
in Fig. 3. It is seen that the order parameter decreases smoothly than in Fig. 5a; 




0.0 0.2 0.4 0.6 0.8 1.0 

2 

a 



Fig. 6. The dependence of the order parameter fj and the susceptibility x vs - noise 
intensity a 2 at e = 0.1, a = 0.1, (3q = 10.0 
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the generalized susceptibility has no pronounced peak. Thus, one can conclude 
that the corresponding phase transition is of the non-critical character. 



6 Conclusions 

In this paper we have considered the physical system of the reaction-diffusion 
kind with a field-dependent mobility and internal fluctuations. It was shown 
that in generalized approach basing on the nonlinear kinetic equation the 
field- dependent mobility is of a power- low form proposed by C.L. Emmott 
and A.J. Bray [11]. Considering the local dynamics of the concentration field, 
it was found that internal fluctuations can sustain stable periodic patterns if 
the noise intensity does not exceed a critical value. The mechanisms of the 
phase transitions have been studied with the help of the mean field theory. It 
was shown that an increase in the noise intensity leads to disordering phase 
transitions. We have found that there is a special domain of the exponent 
determining the power-low form of the field-dependent diffusion coefficient 
when disordering phase transitions are not critical. It can be explained by the 
competition between noise and driving force. The corresponding bifurcation 
and phase diagrams are verified by computer simulations. 
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